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SYNOPSIS 

TEMPERATURE RESPONSE OF WATER BODY WITH 
THERMAL EFFLUENT DISPOSAL 


A Thesis Submitted 

In Partial Fulfilment of the Requirements 
For the Degree of 

MASTER OF TECHNOLOGY 
by 

AJAY KUMAR AGRAWAL 
to the 

Department of Mechanical Engineering 

Indian Institute of Technology/Kanpur 

July 1982 

To predict behavior of prototype systems, the range 
of operating variables for experiments is selected using a 
similitude criteria. Effects of intake location, discharge 
water flow rate and temperature rise are studied for water 
bodies of various depths « A bottom intake is the best 
choice because the buoyancy as well as mixing dominate in 
the upper zone only. Vertical as well as horizontal 
stratification occur during transient regime. At steady 
state an isothermal buoyant pool is established. It's depth 
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is characterized by a non-dimensional pool fraction Which, 
in turn, depends upon the discharge conditions and intake 
location 

Using a heat-balance and corresponding differential 
equations, transient distribution of natural and excess 
temperature in the water body are obtained. Overall heat- 
transfer coefficient at the water surface gives a quick 
guess of the temperature rise at given plant head load and 
weather conditions , Natural water temperature and wind 
velocity are the most influential factors for the determi- 
nation of the excess temperature. 



CHAPTER 


1 


INTRODUCTION 

1 . 1 DESCRIPTION t 

The growing dependency upon electricity for scientific 

developments, industrial growth and prosperity of people has 

called for augmenting the power generation facilities. Todate, 

thermal power plants are the major source of electrical energy 

because of their low initial cost and short gestation period. 

The efficiency of such power plants depends upon the temperatures 

of heat source and sink. While the former is decided by the 

metallurgical limit the latter, in turn, is controlled by the 

choice of waste heat disposal system. Out of about 3 MW of 

energy transfer from the source to the heat engine for a 

typical power plant, nearly 2 MW is rejected to the sink, 

rendering about 1 MW of electrical power output. Therefore, 

water quantity on the order of 40 to 60 litres/s-MW of power 

generation, passes through the condenser, and generally gets 

■fcUe 

heated up by 8-1 2°C, depending upon^ system design. Waste heat 
is then dissipated to the earth's atmosphere through the use 
of cooling systems, that is river, lake, cooling pond, cooling 
tower, etc. 



In once-through cooling systems, water is taken 
from a natural basin, and discharged into another part of 
it, where it is diffused by water currents and subsequently 
dissipated to atmosphere. This situation is favourable only 
for medium and small capacity thermal power plants located 
near sea, river or canal. The condenser cooling water is 
taken upstream and hot water discharged downstream resulting 
in elevated river temperature. It brings about ecological 
disturbances in larger distances. Therefore , the public 
regulatory measures of thermal pollution in river, canal 
etc. may forbid the discharge of hot water from power plants. 

To overcome these limitations, the most common prac- 
tice is the use of closed systems for cooling the condenser 
water. A few available systems are cooling towers, cooling 
ponds, natural reservoirs like lakes etc. In cooling towers, 
heat disposal is accomplisned mainly by ^evaporation of water 
induced by enhanced air-water contact area. In cooling 
ponds/reservoirs (water bodies), the hot water is discharged 
at the surface and cool water withdrawn from the lower level. 
The plant heat is dissipated due to increased evaporation, 
convection and radiation at the air-water interface. 

Because of high initial and operating costs com- 
pared to other closed systems, cooling towers are adopted 
only for zones of expensive land. Hence factors determining 



3 


power plant location invisag.e the utilization of naturally 
available cooling reservoirs till they are exhausted. Such 
systems require minimum operating and initial costs, render- 
ing beneficial results in the long run. Despite involvement 
of additional initial investments, cooling ponds are exploited 
in absence of availability of natural reservoirs. 

Cooling systems employing combinations of once- 
through system, cooling pond/reservoir and cooling tower 
have important influence on optimization among economics of 
power generation, environmental measures and energy conser- 
vation aspects. Adequate knowledge of behaviour of various 
mixed modes of operation require an exhaustive study of indi- 
visual systems. 

The study of cooling water bodies becomes momentous 
in view of their wide spread applicability in single as well 
as mixed mode operation of cooling systems. In the choice 
of a natural lake or design of a cooling pond, a major consi- 
deration is to increase the utilization factor of the basin 
to provide maximum cooling of water and possibility to cool 
water to the steady value of relatively low temperatures. 

In an ideal case, water should reach the natural water 
temperature which, in turn, is governed by meteorological 
and climatological factors of the geographical region, and 
unaffected by warm effluents from the plant. But, in practice 
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a temperature at intake higher than the natural water tem- 
perature is established* This excess temperature/ the 
governing factor of power plant efficiency, can be minimized 
with the knowledge of temperature developments in the water 
bodies subjected to thermal discharge,- under various plant 
load conditions. Cooling efficiency of the system can be 
enhanced by directing the water currents in the cooling 
water body, to ensure maximum utilization of the tempera- 
ture gradients and available surface area. 

To avoid any short-circuiting of warm water, loca- 
tions of intake and discharge points are important. For 
the investigation of behavior of cooling water bodies of 
various depths and configurations, plant variables (that is 
cooling water flow rate, water inlet temperature) and meteo- 
rological variables (that is ambient temperature, relative 
humidity and wind speed etc. ) should be considered. To 
study the ecological aspects of the problem, prediction of 
excess temperature distribution established due to plant 
heat load and meteorological changes is required. 

1.2 LITERATURE REVIEW? 

Different aspects of the problem of heated dis- 
charge into water bodies have been investigated by several 
past researchers. These include the study of heat and mass 
transfer mechanisms at air-water interface, temperature 
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response of water bodies subjected to hot water discharge, 
and jet behavior under various discharge conditions. Expe- 
rimental and analytical techniques have been presented for 
specific eases as reviewed below. 

Raphael [l] developed a procedure including empi- 
rical relations for predicting the temperatures in water 
bodies. The heat and mass transfer coefficients were obtained 
from observed data on lakes. Edinger et al [2] have used 
equilibrium temperature method for the calculation of daily 
average heat-exchange. Since, equilibrium temperature, 
(temperature when net heat flux across the air-water inter- 
face is zero) may vary as much as 40 °C in a day, its use as 
base temperature in linearizing the heat flux equation, is 
doubtful , and may be erroneous . 

Ryan et al [3] developed a technique for account- 
ing the contribution of free convection to the evaporation 
from a cooling pond and putforth its correlations with 
forced evaporation. The variation in the temperature diffe- 
rence encountered in practice is not large. Hence the tem- 
perature dependent term for natural convection, leading to 
complicated expression of heat-flux, can be substituted 
by a constant term. 

Yotsukura et al [4] approximated the basic ex- 
change equations for air-water interface by Taylor series 
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expansion around an undefined base temperature. Error 
analyses were done for linear and quadratic approximations. 

It was suggested that for most practical cases, a linear 
approximation gives acceptable solution. Jobson, H.E.[5] 
used natural water temperature for calculating a surface 
heat transfer coefficient, which, in turn, could be used to 
predict unnatural excess temperature. Kothandaraman, V-[s] 
presented a report upon an investigation regarding nature 
of seasonal and non-seasonal variations in the daily mean 
river temperature at a given location, and using Fourier 
Series analysis developed a method to predict natural water 
temperature on the basis of meteorological data. ' 

Experiments on turbulent water jets discharged 
into a water body were first conducted by Jen et al [?] for 
various jet orientations discharging vertically, horizontally 
and at several angles in between, for different Froude 
numbers. The limitation to the use of results arises as 
the ©tlfice was located nearby at the bottom, which is 
seldom the case. Tamai, et al [8] have determined the mix- 
ing effects of warm water discharge, indicating therein a 
principal role of densimetric Froude number, being played 
in surface spreading of warm water. 

A number of mathematical models, [9, 10, 11, 12] 
right from complicated three-dimensional to simplified 
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isothermal models have been found for prediction of tempe- 
rature distribution in water bodies with space and time. 

Moore et al [ 9 ] have developed a simple model, based on 
perturbation technique, for thermal structure of a strati- 
fied lake. It considers various changes in equilibrium 
temperature and diffusion coefficient to represent heat 
addition and mixing effects of power plants on the thermal 
cycle of the lake. Surface heat exchange coefficient used 
is based on equilibrium temperature method, and its values 
have been taken constant for winter and summer periods. 

Me Guirk, J.J. et al [lo] have described finite difference 
calculations of three-dimensional heated surface jets dis- 
charging into stagnant water. The method presents the solu- 
tion of continuity, lateral and longitudinal momentum and 
energy" equations . Computed results compare well with experi- 
mental data. Similar approaches have been adopted by other 
workers [ll, 12], Most elaborate mathematical studies have 
been presented by Lee, S. et al [l2j. They developed a 
time-dependent, three-dimensional, free surface numerical 
model with horizontal and vertical stretching. The effects 
of variable bottom topography, free surface elevations, 
surface heat transfer currents and meteorological conditions 
were considered. It provides sufficient three-dimensional 
information for study of water bodies. 
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Nystrom et al [13] conducted a number of physical 
model studies , field surveys and analyses. They have empha- 
sized the necessity of similarity analysis between model 
and prototype studies. Froude, Reynold and densimetric 
Froude numbers formed a criteria for similarity analyses. 

Recently the study of cooling ponds has received 
due attention. Jirha et al [ 14 ] presented an analysis of 
thermal stratification of cooling ponds and categorized using 
a non-dimensional pond number, the thermal structure of ponds 
into either the well stratified or partially mixed or ver- 
tically mixed type. Only long term mean thermal structure 
could be predicted. Therefore, any short-term variations 
due to meteorological changes or plant operating conditions 
cause transient deviations from the mean structure. Jirka 
et al [l5] have considered the efficiency of various cooling 
ponds and concluded that deep stratified water bodies are 
most suitable configuration. Otherwise, for shallow cooling 
ponds, a moderate length/width ratio on the order of about 
5-10 as well as some amount of baffling arrangement is 
recommended. 

1.3 PRESENT STUDY ; 

It comprises experimental and analytical study of 
temperature developments in cooling water bodies having hot 
water discharge. Using similitude criteria, the operating 
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variables for experiments are selected with a view to predict 
behavior of prototype systems. Effects of plant operating 
conditions (that is, condenser water flow rate and its tempe- 
rature rise), intake position, pond depth and wind velocity 
on transient and steady-state intake temperatures are studied 
extensively. At steady state, an isothermal layer is estab- 
lished in the upper zone because of turbulent mixing and 
buoyancy. Effects of discharge conditions and intake loca- 
tion on the depth of buoyant pool are investigated. 

An energy balance equation is found for steady-state 
heat-exchange across the water surface. Mass-transfer co- 
efficient, which is a function of wind velocity, is evaluated 
experimentally. A methodology is presented to calculate 
the excess temperature at water surface (at), induced by 
thermal discharge. Meteorological variables, such as solar 
radiation, air-temperature, wind velocity, relative humidity 
etc. influence the thermal structure of water body. There- 
fore, the steady state analysis has been extended to obtain 
the transient variations of A An overall surface heat- 
transfer coefficient, comprising coefficients of convection, 
evaporation and back radiation heat-transfer, is used to 
predict ZA T. Differential equations for natural and excess 
temperature distribution in the water body are obtained and 
solved numerically. 



CHAPTER- 2 


EXPERIMENTAL STUDY 


2.1 THERMAL DISCHARGE INTO WATER BODY; 


Various physical phenomena control the hydrodynamics 
of a heated effluent being transported from the outfall to 
the main receiving water body. In near field region, the 
effluent momentum and buoyancy govern the discharge trajectory 
and dilution. In intermediate region , buoyancy force become 
dominant and control plume spreading and depth. However, in 
the far field region, heat and mass -transfer at the air- water 
interface control the stratification characteristics. 


The relative action of momentum and buoyancy forces 
along jet trajectory is measured by a local densimetric 
Froude number (Fr # ) ; 


where. 


Fr' = V / ( (2.1.1) 

V = horizontal velocity (m/s) 

g' = ( &P / p_^) . g = reduced gravitational 

2 

acceleration (m/s ) 

AP= ( P. - p ) 

. x w 

3 

= density of hot water (kg/m ) 
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p w = density of ambient water (kg/m ) 

Z = depth of water body (m) 

Figure 2-1 exhibits the structure of a buoyant 
surface jet. At near field region-, high values of velocity (V 
and reduced gravitational acceleration render Fr' >> 1, since 
( Ap / p_^) ± s nearly between 0.002 - 0.005 giving g'Z << 1. 
Momentum causes linear growth of jet, both laterally and 
vertically, and entrainment of surrounding fluid into turbu- 
lent jet region. 

At larger distances Fr' decreases as V and g' decrease 
asymtotically. Therefore buoyancy exerts an increasing in- 
fluence on the jet behavior, and starts dominating the inertial 
forces. Finally it causes an unsteady motion in all direc- 
tions, rendering a gradual increase in length, width and thick- 
ness of the pool of buoyant water. After a long period of tim< 
a steady state, based on a balance between heat-loss to the 
atmosphere and heat addition through the discharge, results 
in the establishment of a buoyant pool in all the regions. 

2.2 SIMILITUDE CRITERIA: 

The study of prototype is, in general, extremely 
complex. Therefore it is a usual practice to simulate the 
prototype system using a model . The governing dimensionless 
numbers are sought to establish similarity. In the present 
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type of problem, model and prototype Froude number (Fr) 

.should be equal. It represents the ratio of inertia force 
to gravity force, that is, 

Fr = V / ( v/g Z ') (2.2.1) 

for similitude 

Fr = Fr (2.2.2) 

m p 

where subscripts 'm' and 'p' are used for model and proto- 
type, respectively. 

From Eqs. (2.2.1 and 2.2.2) with respective quantitie 
one gets ; 

V r = /Z ' (2.2.3) 

having V = V /V and Z = Z / Z 
^ r nr p r m p 

For kinematic similarity Reynold number 'Re' 

( = V Z / v , where v is kinematic viscosity) which measure 

the relative action of inertia to viscous force, is used. It 

governs the flow pattern and mixing process. Therefore, 

Re r " Re n> / Re p = V r Z r / v r (2 - 2 - 4) 

Since v is unity, using Eq. (2.2.3), Re m can be 
expressed as s 

Re = Re (Z ) 3//2 (2.2.5) 

m p r 
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Prototype Reynolds numbers are well within turbulent 

regime. Therefore, the value of model Reynolds number depends 

upon depth ratio . A minimum limit to the value of Re_ is 

m 

imposed (500-1000 ) such that the model operates in turbulent 
regime [ 15 ] * This subsequently sets a lower limit on Z^_ . 

2T 

To achieve similitude of buoyancy, the densimetric 
Eroude number must be equal in model and prototype. Therefore , 

Fr r = v r 7 ( / "Up /P ± T r 7g r ”J z r "0' = i ( 2 . 2 . 6 ) 

Using Eg. (2,2.3) into Eq* (2.2.6) with g = 1, 
one gets ; 

( A P / P ± ) r = 1 (2.2.7) 

Eq. (2.2.7) can be satisfied only if temperatures in model 
and prototype are the same. It implies , therefore, that for 
equal temperature rise at inlet, temperature profiles for 
model and prototype become similar. Therefore, the model 
study would enable the prediction of temperatures in prototype 
systems . 

2 . 3 EXPERIMENTAL SET-UP 

It consists of a water tank with various measuring 
and control devices for investigation of temperature deve- 
lopments and prediction of important design parameters. 

The arrangements for the study of effects of wind velocity 
and location of inlet and discharge points have also been 
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made. The overall set-up encompasses broader flexibility 
to facilitate the experimental study at desired conditions. 

A schematic diagram of the system is shown in 
Figure ( 2 . 3 ) and photographic view in Fig. (2.2(a) to (c)). 

A componentwise description of the experimental set-up is as 
follows s 

2.3.1 WATER STORAGE TANK 

The water body is simulated by a tank of size 570 x 
570 x 425 mm deep, fabricated from 12.7 mm thick perspex sheet 
For discharge and intake, holes were drilled at various posi- 
tions as shown in Fig. 2.4(a). Intakes not under use were 
kept closed. The tank walls and bottom were insulated with 
25 mm thick thermocole sheets with a view to minimize heat 
losses from walls and bottom to surroundings. 

2.3.2 WATER CIRCULATION ARRANGEMENT 

Warm water was discharged at the top of the water 
body having cool water intake at the lower level. Teflon 
tubes, insulated with thermocole, were used for all the connec- 
tions. The steady flow of hot water in the discharge line was 
achieved by an adjustable constant head device. Fig* 2.4(b). 
The flow rate was measured with the help of a calibrated 
rotameter. It is worthwhile to mention that the direct 
control of steady hot water flow rate, even with the help 
of most sensitive valve was found to be unreliable. At 
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intake, the rotameter was mounted in between intake point and 
constant temperature bath (kept at level lower than the 
tank) , The desired flow rate, equal to discharge rate was 
achieved by adjusting a Hoke valve. 

To simulate with hot. water circuit of thermal power 
plant, an ultrathermostatic bath (Temperature range: - 30 to 
200°C ) capacity s 15 litres, control accuracy : + 0.03°C ) 
heating capacity ? 1500 W) was used. The desired temperature 
of the warm water could be achieved by adjusting a mercury-in- 
contact thermometer having an on-off heater control. , 

2.3.3 ARRANGEMENTS FOR TEMPERATURE MEASUREMENT z 

Temperatures were recorded using 20 gage copper - 
constanton and 36 gage Alumel - chromel thermocouple wires. 
Ordinary copper wire was used for extension from thermocouples 
to calibrated temperature indicator potentiometer. 

For the study of stratification in the water body, 

36 thermocouples were fixed in four horizontal planes, each 
plane having 9 thermocouples. In addition, 15 thermocouples 
were placed uniformly at the water surface. Fig. 2.4(c). 
Separate -thermocouples were used to get intake, discharge 
and atmospheric temperatures. 

A moving probe, exhibited as in Fig. 2.4(d) consists 
10 thermocouples glued to two 2 mm thick copper sticks. 
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It was mounted on a sliding stand having movement along the 
jet direction in the horizontal plane. The probe itself 
slides laterally (perpendicular to jet direction) and verti- 
cally on this stand. This arrangement facilitated the mea- 
surement of temperature at any desired location in the tank* 
Measuring scales were fitted to locate any point of interest 
in the tank. 

2.3.4 TEMPERATURE INDICATOR POTENTIOMETER 

A 48-point Honeywell temperature indicator poten- 
tiometer was used for successive temperature measurements. 

It was pre-set for copper- constanton thermocouples and had 
a measuring accuracy of 0.15 °C. The indicator scale was, 
however/ calibrated for alumel-chromel thermocouples. 

A 20-point temperature recorder was used to record 
continuous variations of surface temperature at the beginning 
of the experiment. The values were recorded graphically 
within an accuracy of + 1°C. 

2.3.5 SURFACE AIR DISCHARGE ARRANGEMENT 

It consists of a blower mounted at the upper level 
of the tank. The air was supplied through a round duct 
having ractangular exit (600 x 130 mm deep) placed near the 
tank. To have uniform air stream, wire-mesh strainer was 
put in the duct. The desired wind velocity (in the range of 
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4 to 16 Km/h) was maintained by regulating the supply voltage 
of the blower » 

2,4* EXPERIMENTAL PROCEDURE; 

2.4.1. RANGE SELECTION 

To simulate water bodies having usual depth 15 to 
30 m , depth ratio (2 ) is taken to be 1/30 to 1/60, Subse- 
quently velocity ratio should vary between 0.1 3 to 0*18* 

For prototype systems * Reynolds number (Re ) is of the order 
6 

of lo and discharge velocity vary between 0.2 to 0.4 m/s. 

The ranges for model Reynolds number (Re m ) and jet velocity 
are found to be 2000 to 6000 and 1.1 to 3.0 m/s, respectively. 
Hence, jet diameter (d ) is taken to be = 1*3 mm and hot 
water flow rate 100 to 200 ml/min. 

Densimetric Froude number limit the temperature 
rise at intake ( Therefore, bT^ is selected as 10 to 

20°C, the normal range of condenser temperature rise in a 
power plant, 

2,4,2 PROCEDURE 

First, the storage tank is filled upto desired 
level with water at ambient temperature. Heater in ultra- 
thermo-static bath is put on, and murcury-in-contact-thermo- 
' meter adjusted to give preselected temperature rise at the 
discharge point. 
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Height of the constant head device is adjusted to 
yield desired hot water flow rate, as checked by rotameter. 

The same flow rate is maintained at intake by proper opening 
of Hoke Valve. 

All the temperatures, inside the tank and at the 
surface of the water, are measured and recorded at 1/2 to 
2 hour interval, depending upon the case. Measurements are 
continued till a steady-state, which occurred in 6 to 12 hours, 
is achieved. 

Variations in intake temperature and depth of 
buoyant pool are observed for a number of discharge condi- 
tions and intake locations. Finally, experiments are con- 
ducted for a set of wind velocities to determine heat and 
mass-transfer coefficients. 




CHAPTER- 3 


ANALYTICAL STUDY 

3.1 DESCRIPTION: 

At steady state meteorological conditions, the 

water surface eventually reaches an equilibrium temperature 

(T ) for which the net heat-flux at the interface is zero, 
e 

But in practice, the water body is subjected to rapid varia- 
tions in meteorological data. Hence its equilibrium tempe- 
rature has large fluctuations even for a single day. How- 
ever# it always try to reach T asymtotically with time due 
to thermal inertia. Therefore# the actual temperature 
attained by the water body# called as natural water tempe- 
rature (T )# has very small fluctuation compared to T . 

n *== 

The temperature distribution of the water body 
alters in presence of a heat source. This gives rise to 
an excess temperature (^1) over T n # and is determined with 
the help of energy balance equation. 

3.2 ENERGY BALANCE EQUATION 

The net short-term heat exchange across a unit ■ 
area of the surface per unit time# Q (T#t)# can be expressed 


as : 
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Q (T,t) = Q q (t) + Q a (t) - Q b (T) 

- Q e (T,t) - Q c (T,t) (3.2.1) 

where 

T = surface temperature of water (°K) 

Q G (t) = heat gain by short-wave radiations 

2 

absorbed by water at time t (W/m ) 

Q a (t) = heat gain by long -wave atmospheric 

radiations absorbed by water at time t 

(W/m^ ) 

Qb(t) = heat loss due to back radiations at 
time t (W/m^ ) 

Q (T,t)= heat loss due to evaporation (W/m ) 

W 

2 

Q c (T,t)= heat loss due to convection (W/m ). 

Q G (t) is directly known from radiation data. 
Remaining quantities can be obtained in terms of natural 
water temperature and meteorological variables as follows: 

3.2.1 ATMOSPHERIC RADIATION 

The estimation of atmospheric radiation is extre- 
mely complex due to involvement of moisture/ temperature, 
ozone, carbon dioxide etc. However, using weather obser- 
vations, it is expressed as [l]s 

Q a (t) = 0.97- a. J3 . T a (3.2.2) 

where 0 « 0.75 to 0.85 
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3.2*2. BACK RADIATION 


Long wave radiations emitted by a water surface 
can be computed from the well-known Stef an-Boltzman law 
for the black body radiation : 

Q b (t) = 0.97 . CT. t 4 (3.2.3) 

3.2.3. EVAPORATION 


It is caused due to the difference between the 
vapor pressures at water surface and ambient air. Its 
magnitude depends upon air temperature, humidity and wind 
speed etc. The heat- transfer by evaporation is, however, 
evaluated from 


C3 (T,t) = h(W )* [e(T ) - E (t)] . L(T) (3.2.4) 

0 3 


where , 

h(W) 

W 

E(T) 

E (t ) 

3 

0 

L(T) 


= mass-transfer coefficient = (a+bW) 
'a' and 'b' are constants 
=s wind velocity (km/h) 

= saturated vapor pressure of air at 
temper a tur e 'T' ( bar ) 

= vapor pressure of ambient air (bar) 
= 0 , E(T ) 

cL 

= relative humidity 

sa enthalpy of evaporation at tempera- 
ture 'T' . (KJ/kg ) 
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3.2.4. CONVECTION 

On the basis of heat and mass transfer similarity 
the ratio Q / Q has been given by Boven [l6] as s 

C 0 

Q / 0 = 0.61 x 10" 3 -(T-T ) / [E(T) - E ] (3.2.5) 

C ' "G a. cl 

hence, 

Q (T,t) = 0.61 X 10"' 3 . h(W) . [T-T (t)] . L(T) 

C a- 

.... (3.2.6) 

The overall expression for heat-exchange at water 
surface takes the forms 

Q (T, t ) = Q Q (t) + 0.776 . a . [T a (t)] 4 - 0.97 . a . T 4 

- h(W ) * [E (T ) - 0 . E(T a )] . L(T) 

- 0.61 x 10~ 3 .h(w) .(T-T ) . L(T) (3.2.7) 

a 

with, 

E(T) « .02 33.8 EXP [l8.1 - 5303. 3 /t] (3.2.8) 

L (T ) = (3145.99 - 2.368 . T) . 10 3 (3.2.9) 

3.3 EXCESS TEMPERATURE (AT): 

Heat balance equation is separated into two compo- 
nents that is, first term represents heat— flux caused natu- 
rally while the other accounts for heat loading. Using 

T =; T + at, one can express: 
n 

Q (T,t) = Q (T n + AT,t) 

T +AT 

= Q (T , t) + Q t (T,t).dt (3.3.1) 

n 
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where , 

Q (T , t) = heat-transfer at natural water tempera- 
ture , T 

n 

Q T (T,t) = £ Q (T,t)/'5T (3.3.2) 

Second term on R.H.S. of Eq, (3.3.1) may be inte- 
grated numerically using Trapezoidal rule as % 


T 4AT 

Q t (T,t) . dt = . [Q t (T n , 


For a preliminary estimate of aT, 
linearized taking T n as the base temperature 


t) + Q t (T n +AT,t) 

.... (3.3.3) 

Eq. (3.2.7) is 


Q (T n + AT, t) = Q (T n ,t) +AT . Q t (T n , t) 

.... (3.3.4) 

Using Eq. (3.2.7), Eq. (3.3.4) is expressed as 

Q (T n +AT,t) =s Q (T n , t) - 3.88 . a . T 3 

- h(W) .[L(T n ).E T (T n ) + 

L t (T n )„[E (T n ) - 0 . E (T a )]].AT 

- 0.61 x 10“ 3 . h(W X [l (T n ) 

+ L t (T n ).[ T - T a ( t ) ] ] . AT (3.3.5) 

where , 

E (T ) = ^E(T)/^t/ = 5303.3 . E(T)/T 2 

x u n 

«»•« (3*3.6) 


# 



24 


LT ( Tn )= fe L tI )/ )T f T=Tn = . 236S (3 . 3 . 7) 

Using an overall surface heat transfer coefficient 

2 

(K; W/m -K ) , Eq. (3.3.5) is expressed as 


Q (T n + AT, t ) = Q (T ,t) - K . AT 


(3.3.8) 


with. 


K = K (T n ,t)= - Q t (T ,t) = + K e + k c (3.3.9) 

= back radiation heat-transfer coefficient 
,3 


= 3.88 . a . T. 


n 


(3.3.10) 


K 0 = evaporation heat-transfer coefficient 

= h(W) [L(T n ).E T (T n )+L T (T n X[E(T n )-0.E(T a )]] 


(3.3.11 ) 


K c ss convection heat- transfer coefficient 

= 0.61 x 10“ 3 . h(W). [l (T )+ L_ (T ) j*T— T ( t ) ] ] 

... (3.3.12) 

Equation (3.3.8) gives the value of excess temperature as: 

AT = QP/K (3.3.13) 

where , 

QP * plant heat load (W/m 2 ) = Q(T n# t) - Q(T n +^T,t) 

Excess temperature predicted by Eq. (3.3,13) is 
used together with Eq. (3.3,3) to iterate for actual tem- 
perature rise ( at' ) 

AT' = 2 . QP/[K(T n ,t) + K(T n + AT, t ) 


(3.3.14) 
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At known natural water temperature and meteorological con- 
ditions , Eqs. (3.3.13 and 14) are used to determine tempe- 
rature rise at the surface for the given plant heat load. 


3.4 DIFFERENTIAL EQUATIONS AND METHOD OF SOLUTION: 

In developing the differential equations following 
assumptions are mades- 

1. Energy transfer across sides and bottom of the 
water body is insignificant. 

2 . Temperature variations occur only in the vertical 
direction* Horizontal gradients are negligible. 
They are adequately supported by the experimental 
results . 


3. Water body is rectangular. 

4. Properties such as p , c» etc. are constant. 

IT 

3.4.1 DIFFERENTIAL EQUATIONS: 

The thermal-energy stored by water, relative to 
that contained at a base temperature (T 0 ), is given: 

Energy .= c^ , P [ (T n + at) - T 0 ] 

The conservation of thermal-energy equation for 
a differential volume of water is ] 

^[(T n + AT) - Tg] V[ ( T n + aT) - T b^ 


(3.4.1) 
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that is. 


^ T n T B . 'b(^T) ~d ^ T n T b\ ^ 2 (AT) 

+ ~ = s 7 — r — 5 + ‘ 

z d Z 

(3.4.2) 


dt 


<*t 


Boundary conditions: 


At 


At 


Z = 0 : e 


^(T n 4AT T b ) Q(T,t) 


"6 Z 


P • c r 


^(T + AT - T ) 

Z=H : £ Z— ^ - = ° 


(3.4.3) 

(3.4.4) 


An equation similar to Eq. (3.4.1) is derived 
for the distribution of natural water temperatures : 


”&t 

Boundary conditions 
At Z = 0 ? i 


(T n ' V 

~ = £z ~T? 


<>«„ - V Q (T n , t) 


'z 


T> z 


P • c 


With Q (T ,t) from Eq. (3.2.7) and 


At 


Z = H 


(T - T ) 
v n 3 


= 0 


(3.4.5 ) 


(3.4.6) 


(3.4.7) 


Since the natural water temperature distribution 
is independent of the excess temperature/ the distribution 
of at can be had by substracting Eq. (3.4.5) from Eq. 
(3.4.2 ) . It gives : 


AT) 

*"b t 


V < AT> 

" 

O Z 


(3.4.8) 
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With Boundary Conditions: 


At 


and 


Z = 0 


Z = H 


£ z [ O ( AT) / b z] = K 
AT 


AT /Pc 

P 

... (3.4.5) 


e 2 [ b>< 4T) /bz ] = o 


(3.4.10) 


3.4.2 METHOD OP SOLUTION: 


Equations (3.4.5 and 3*4.8) are solved numerically 
to get natural and excess temperature distribution. Heat- 
transfer across the surface is calculated using hourly 
meteorological data. The cooling water body is divided 
into 'n' horizontal elements, each of depth AH. Finite 
difference equations are written for each element. 

(Appendix-I ) . Since time step is one hour, the necessary 
requirement of convergence is fulfilled by a proper selec- 
tion of depth step. Eq. 3.4.5 is solved by iterative tech- 
nique for hourly natural water temperature distribution in 
the water body. Iteration is continued until the conver- 
gence is obtained, indicated by a small difference (0.01 to 

‘P+l 

0.1°C) between T n (I) an d t through out the year, for 

all values of I considered, (p is the number of iteration 
and I is the number of the depth element from the top). 
Initial guess of T^ (i) is found to affect convergence 
significantly. If the water body is taken to be at cons- 
tant temperature and boundary conditions are imposed 
suddenly, there will be high values of flux near the surface. 
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Therefore, the values of ^ T n /<> z computed during the 
first few steps in time will be very inaccurate. Hence 
instead of constant values at t = 0, the iteration -is started 
from a computed solution at t = 10 to 12-hours. 

Excess temperatures at the surface throughout the 
year are calculated using T n at the surface and Eqs. (3.3.13 
and 14). Eq. (3.4.8) is then solved by iterative technique. 

A linear excess temperature distribution as an initial 
guess is found to give rapid convergence. 



CHAPTER -4 


RESULTS AND DISCUSSION 


Experimental investigation has been carried out to 
study the temperature characteristics at intake with respect 
to its location, pond depth, plant and meteorological vari- 
ables for transient and steady conditions. Following dimen- 
sionless quantities are used to represent variables; 


Dimensionless 


intake temperature; T q = (T 


2 

Dimensionless time r, ^ = or, t / H 
where , 


T q = intake temperature ( °K ) 

T a = atmospheric temperature (°K) 
T^ = inlet temperature (°K) 
t = time ( s ) 

H = pond depth (m) 

2 

a = thermal diffusivity (m /s ) 
4.1 INTAKE TEMPERATURE ( T ) 


(4.1) 

(4.2) 


4.1.1. EFFECT OF INTAKE POSITION 


Figures 4.1 and 2 indicate the transient behavior 
of temperature at intake for its various positions . It (T ) 
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rises exponentially with time, having lower value at deeper 

intake and vice-versa. The difference in T values of two 

o 

consecutive intakes in a vertical plane is observed to de- 
crease with decreasing depth. Therefore, a deeper intake 
should be preferred. 

At higher flow rate (or Re), the temperature rise 
is more at all levels. The effect becomes pronounced towards tl 
free surface. It is because the mixing in the upper region 
is much more at higher Re. The experimental results suggest 
that the mixing is, further, accelerated by placing intake 
in this region. 

4.1.2 EFFECT OF HOT WATER FLOW RATE (m) 

The variation of intake temperature, at different 
flow rates and discharge temperature rise a ay = 10 and 15 °C 
for a typical intake position, is exhibited in Figs. 4*3 and 
4. It is observed that intake temperature augments with 
increasing ifi. The change in steady temperature with an 
increase in flow rate from m to 1.5 m is double than that 
for increase in m from 1.5 m to 2.0 m . It is because of 
more mixing at higher flow rate. 

4.1.3 EFFECT OF DISCHARGE TEMPERATURE RISE (AT^ ) 

Transient variations of intake temperature for a 
flow having Reynold number of 4900, and = 10, 15 and 20°C 
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are shown in Fig. 4-5. The increase in aq\ does not pro- 
duce a proportionate increase in T q . It is because the heat 
dissipation across the water surface increases due to higher 
A T^ , 10 to 15 °C, resulting in a corresponding low tempe- 
rature rise at intake. On further increase in at^ , that is, 
from 15 °C to 20°C, the available surface of the water body 
is not adequate to dissipate the heat proportionately. It 
results in elevated intake temperature due to heat accumula- 
tion. Therefore, AT\ more than 15 °C for such a water body 
( corresponding to a heat loading of 600 T - 7 /m or more) should 
not be employed. However, for the same heat discharge, a 
temperature rise between 10-15°C and correspondingly higher 
flow rate is a better choice. 

Figure 4*6 exhibits an overall summary of steady 
state intake temperatures under various discharge conditions. 
It is evident that a higher temperature rise (between 10 - 
15 °C) and Reynolds number (Re m = 3000 to 5000) are a suitable 
choice if regulatory measures of temperature rise are within 
satisfactory limits. For Re m less than 3000, at^ should 
be kept below 10°C. 

4.1.4 EFFECT OF POND DEPTH 

In practice, several types of water bodies having 
different shapes and depths are encountered. However, 
depth is an important parameter. Fig. 4»7 shows the effect 
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of pond depth on intake temperature. It is seen that in 
shallow ponds temperature increases quickly and eventually 
becomes almost uniform. It implies, therefore, that intake 
can be located anywhere except in the vicinity of the dis- 
charge point. The above results have been found to be true 
for all the operating ranges. 

4.2 POOL FRACTION ( ip ) 

Experimental results reveal that a buoyant pool .'is 
established in the upper part of the water body. The 
temperatures below it vary mildly and seldom exceed 10/ of 
intake temperature rise. For practical purposes, a water 
body can be divided into two distinct layers, the upper one 
being warmer and lower one cool. The variation of pool 
fraction (IP : expressed. as the ratio -of warm layer to the 
total pond depth) for various discharge conditions and 
intake locations are displayed by Fig. 4*8[l7] • Depth of warm 
layer increases' with increasing flow rate (or Re ) due to 
enhanced momentum causing mixing. The position of intake 
point has important influence on the value of IP. At bottom 
intake, IP decreases with increasing A T^ while at top intake, 
the behavior is reverse. Temperature for bottom intake 
is comparatively low due to longer flow path. Also the 
higher value of , A T^ causes larger heat dissipation from 

the water surface. It results in proportionately low tempe- 

* 

rature build-up of the water body, rendering a reduction 
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in IP. At top intake, T q which corresponds to the surface 
temperature, is quite high. A major portion of heat supplied 
is carried by the intake water. It decreases the heat loss 
from the surface. The water body is thus subjected to heat 
accumulation which, in turn, results in deeper warm water 
layer. This tendency becomes pronounced as the temperature 
rise at discharge increases. 

4.3 EFFECT OF WIND VELOCITY (W ) 

The steady state structure established in a water 
body gets disturbed due to meteorological variables .like 
wind velocity, air temperature and relative humidity etc. 

A typical nature of intake temperature variation with wind 
velocity is exhibited in Fig. 4»9 for Re m = 5000. It shows 
exponential decreament in T q with increasing W. The graph 
signifies the importance of considering meteorological 
parameters for the analysis of water bodies. 

Since experimentation is not possible in transient 
environments, an analytical solution is required. At zero 
wind velocity, the natural convection and diffusion govern the 
surface temperature. Heat- transfer across the interface 
increases at higher wind velocities . Surface temperatures 
are obtained experimentally for various cases. Heat and 
mass transfer coefficients are then calculated from energy 
balance equation. Fig. 4.10 shows the variation of mass 
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transfer coefficient with wind velocity. A linear profile 
is selected. Values reported by Yotsukura et al [ 4 ] are 
low as effect of natural evaporation has been neglected. 
h(W) obtained by experimental results match well with its 
values reported by Ryan et al [ 3 ]. 

4.4 ANALYTICAL RESULTS 

Using hourly values of meteorological data of 
N. Delhi [l8], transient variations of daily average natural 
water temperature are presented in Fig. 4*11. The changes 
in temperature with depth are not large because of mixing 
caused by wind currents etc. 

In Fig. 4*12, transient variations in excess tem- 
perature at the surface (at) are shown. at is minimum 
in June, July and maximum in December, January. It is because 
during winter, T n is low. At low temperature, heat dissipa- 
tion from the surface is less. Therefore, to dissipate the 
plant heat, the surface temperature should elevate. Reverse 
is the case for summer. As the heat load increases, at also 
increases, but not proportionately> as seen from Fig. 4*12. 
Error of linearization is more for higher AT but for a 
temperature rise of few degrees (2 to 4°C) it is almost 
negligible. 
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Figure 4 » 1 3 exhibits the transient distribution of 

excess temperature in a water body, for a plant load of 
2 

lOO W/m . Sharp gradients occur near the surface and compa- 
ratively less fluctuations are observed at lower depths. The 
net temperature distribution can be had by summing up T n and 
AT due to a certain plant load, at each location. 

Excess temperatures are predicted using an overall 
surface heat- transfer coefficient (K) defined by Eq. (3.3.13). 
Fig. 4«l4 shows the effect of meteorological variables and 
T n on K. It's values for two limiting cases of air-tempera- 
ture and relative humidity (T = 50°C, 0 = 1000 and T = 0°C, 

3 3 . 

0 = 0/) are presented. It is evident that T and 0 have 

3 

almost no influence on the value of overall surface heat- 

transfer coefficient. Therefore, K may be expressed in terms 

of wind velocity and T^ only. Generally T^ vary from 10 to 

40 °C and W from 6 to 20 km/h. Hence, a practical zone of K 

2 

values (varying between 15 to 90 W/m ~K) is identified. 

K consists of back radiation, evaporation and con- 
vection heat- transfer coefficients. It is exhibited from 
Fig. 4*15 that evaporation is the dominating cause of heat- 
dissipation. Moreover, it' changes largely with wind velocity 
and T . 
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4.5 ILLUSTRATION” 

In the following illustration, present results 
are used to select the cooling water system of a thermal 
power plant . Its effect on temperature rise in the water 
body is also obtained. 

Available Information:: 

1 . Ultimate generation capacity of the thermal power plant 
complex s 5130 MW* using 210 and 500 MW steam turbines. 

2. Condenser water cooling : Through the reservoir of 

a hydel power plant located 
nearby. 

Reservoir Data; 

1 . Surface Area : 


(a ) 

At 

maximum 

water 

level, s 

70, ,000 acre or 

2.83 x 

1o 8 2 
10 m 

(b) 

At 

minimum 

water 

level ; 

32,000 acre or 

1.29 x 

o 

00 

3 

to 


2 . Water Storage s 

(a) Live storage at max. water level ; 7.28 x 10 acr— ft 

or q - 
8.98 x 10 rri 

6 9 3 

(b) Dead storage; 1.32 x 10 Acr~ft or 1.63 x 10 m 
SOLUTION; 

Rate of heat conveyed by cooling water; 2 x 5130 

= 10,260 MW 

* 

Therefore, m = 2.45 x 10 3 kg/sec.for A T^ = 10°C 
= 1.64 x 10 3 kg/sec.for = 

Plant heat load l QP = 35 to 80 W/m 3 


15°C 
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This gives, 

2 = 58.8, V = 7.67 

r r 

and 

Re m = 2.216 x 10~ J . Re 
m p 

Since Re^ = 3000 to 5000 with a temperature rise of 

10-15 °C is preferred, the desirable Re is obtained ass 

P 

Re p = 1.35 to 2.25 x 10 6 

let diameter of each discharge pipe = 5 m. 

Therefore, discharge velocity = 0.24 to 0.4 m/s 

discharge flow rate= 4.77 to 8.0 x 10 kg/s 
Since heat supplied per unit area of the water body is low, 
the temperature rise of 15 °C with a lower water flow rate 
shall be a suitable choice. Discharge channels can be placed 
correspondingly. 

Now, 

Natural water temperature s 17.5 to 31 °C 
(from Pig . 4.il ) 

Wind velocity s 6 to 13 km/h 

Therefore, overall surface 

2 

heat- transfer coefficient : 18 to 60 W/m -k 

Approximate excess temperature: 0.58 to 4.4°C 

Actual excess temperature 0.5 to 4.0°C 

(using Fig. 4.12) 

Temperature rise shall be maximum during winter. 
Since, excess temperature is well within acceptable limits, 
ecological problems shall not arise. 



CHAPTER-5 


CONCLUSIONS 


From the present study carried out experimentally 
and analytically the following conclusions are arrived at; 

1. A deeper intake allows accelerated mixing at top 
layer/ avoids short circuiting and gi Ves better 
intake tempera tur e . 

2. The higher temperature rise 10 to l5°c anc j p e 

m 

(3000 to 5000) permits effective Utilization 

of available surface area of the wate r body having 

,2 

thermal loading of 200 to 600 W/m . 

3. In shallow water body, almost uniform water tem- 
perature is attained. Therefore/ the position of 
intake is not very critical except i n the vicinity 
of the discharge point. 

4. A steady state buoyant pool is established depend- 
ing upon the discharge conditions, tempera*- 

tures in the upper zone are almost constant, while 
the thermal structure below it remai ns nearly 


unaffected . 
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5 . The wind velocity causes sharp decrease in the 
intake temperature. Therefore, it is important 
to consider its effect on T^. 

6. Heat and mass -transfer coefficient depends linearly 
upon wind velocity, 

7. The natural water temperature in a water body can be 
calculated ' >f rom the meteorological data. 

8 . The excess temperature can be predicted toy lineariz- 
ing the energy balance equation with natural water 
temperature as base temperature. Error of linear! ~ 
zatlon is negligible ^temperature rise of 4-5 °C« 

9. Overall surface heat- transfer coefficient, which 
depends mainly upon the natural water temperature 
and wind velocity, gives a quick guess of excess 
temperature at the water surface. 
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FIG- 4 6 VARIATION OF INTAKE TEMPERATURE 
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FIG-4-8 VARIATION OF POOL FRACTION WITH 
DISCHARGE CONDITIONS. 
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APPENDIX- I 


Eg. 3.4.5 is written as 


?;t 


^b 2 T 

■£z 2 


(A . 1 ) 


where 


T = T - T_ 
n B 


12! _ _2. ( T — T ) 

S>t At ■ i,j+l i , j 


(A. 2) 


_ j_ , 

’ O ~ AT T v - 1 4 

>Z 


^ - -ttt v j. . . - . - 2T . . + T. . .) 
2 Ai l+l , ] i,j l-l, J 


(A. 3) 


This gives 


T 

i , j +1 


At . e r 


ah 


[ T. , . . + T. - .] 

L l+l, J l-l, J 


Cl - 


2a t .e 


Z ] T 


AH 


1/ J 


Using M = 


ah 


2 At . e r 


(A. 4) 
(A. 5 ) 


T, .... = M [t. , . + T . ] + [l - M] T. . 

l,j+l L 1+1, J l-l, J i/J 


(A. 6) 
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Stability condition for A. 6 require : M < 1/2 
Finite difference representation of A. 6 is': 


T 

T 


1 , 

2 , 


J+l - M ^ T 2,J + T 0,J^ + 


J+l = M (T 3,J + T 1,J } + 


(1 - 2M) Tj j 
(1- - 2M) t 2#j 


A.t 


T 

T 


n, J+l 


= M (T 


n+1 / J+l 


= M 


+ T 

n+1, J n-1, 

(T , „ T + T 
n+2 # J n,J 


j> + <! - 2») T n,J 
) + (1 - 2K) T n+liJ 


Temperatures at fictitious points, T Q j and T^ +2 j 
are obtained from the finite diffeJ3pp&e* equations of the 
boundary conditions (Eqs. 3.4.6 and 7) which are: 

(T 2' J+ 1 - T 0.J + 1> f 2 ^ ■ ‘ Q (T 1,J' At, / f c p £ Z 


and T n+ 2 # j “ T n ,J “ 0 


A. 8 
A. 9 


Therefore, the finite difference representation of various 
elements in the water body is: 


For 


I = 1 : T t = 2m[t i+1 ^ j + ah Q(T t t , J .At)/P C^ ej 


I, J+l 


I,,J' 


'p Z" 


+ (1 ~ 2M) T t _ ••• A. 10 

1 < I < n+1 : Tj ^ J+1 = M [ T I+ i # j + T i-i, j^ 

+ (1 + 2M) T_ T ‘ A.1T 

J- !► w 

1 « n+1 : = T i-l,J + ^ T I, J A * 12 

Similar difference equations are written for eq. (3.4.8) also. 



